Crystallization in single- and multicomponent neutron star crusts

We use molecular dynamics simulations to study the formation and stability of single- multicomponent lattices in outer crust neutron stars. Including an improved treatment for Gaussian charge distributions ions, we obtain expressions potential forces arising from electron screened Coulomb interactions using efficient Ewald sum procedure. Our findings show that. baryon densities crust, a pointlike ion cannot fully describe crystallization behavior. In our work, usual parameter, ${\mathrm{\ensuremath{\Gamma}}}_{C}$, along with screening parameter polarizable background, $\ensuremath{\kappa}$, are complemented additional $\ensuremath{\eta}$, providing information on finite size ions. approach find that including beyond-pointlike approaches plasmas under Thomas-Fermi approximation has strong impact calculated lattice energetic stability, decreasing energies per up $\ensuremath{\approx}40%$ respect interaction and, as consequence, melting point is displaced lower temperatures.